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Abstract 

Using a recent global analysis result after the precise measurement of 6n, a possible Herimtian neutrino mixing 
ansatz is discussed, the mixing matrix is symmetric and also symmetric with respect with the second diagonal line in 
the leading order. This leading order ansatz predicts #13 = 12.2°. Next, consider the hierarchy structure of the lepton 
r*j ■ mass matrix as the origin of perturbation of the mixing matrix, we find that this ansatz with perturbation can fit current 



data very well. 
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I. INTRODUCTION 



Many neutrino experiments have revealed that the three known light neutrinos must have finite but small 
masses and that different flavor neutrinos oscillate from one to another. The massive neutrino states vl are 
relateded to the flavor states v/z, by the 3 x 3 PMNS matrix: 



v, 



UpMNS 



I \ 

V3 



(1) 



>L v ' L 

The PMNS matrix is unitary, and can be parameterized in the standard way, with three mixing angles and 
three CP-related phases: 



Upmns - UK 



(2) 



with 



U = 



C12C13 



S12C13 



s l3 e 



-C23S12 - S2}C\2Sne lS C23C12 - ^23*12*13^ ^23^13 

^ 523*12 -czicnsne' 6 -szscn- cnsnsvse® C23C13 



K = diag(e ip ,e icr , 1), 



(3) 



(4) 



Where c;y = cos Oij and Sjj - sin 6^, the phases p and cr are physically relevant only if the massive neutrinos 
are Majorana particles. 



The mixing angles 623 and #12 are much larger than #13. Before the recent experiment results of T2K [1], 
MINOS [2] and Double Chooz [3], the best-fit of neutrino oscillation angles are usually taken to be #13 « 0°, 
#23 * 45° and On ~ 34°. Many ansatzes are proposed to explain the smallness of #13 and maximality of 
sin #23 [6]. The well-known tri-bimaximal mixing pattern is consistent with experiment data then [7]. But 
Daya Bay SBL experiment made a precise measurement of #13 from the v e — > v e oscillations [4]. The #13 
best-fit (±lcr range) result is 

sin 2 29 l3 - 0.092 ± 0.016(stat) + 0.005(syst) , (5) 

which is to say that 9 U t 0° at the 5.2cr level. The RENO collaboration confirmed the Daya Bay result 
soon [5]. 



013 is not really very small. This fact makes many ansatzes face difficulties [6] and opens the door of 
possible new phenomenological applications. The possible explanation and impact of large #13 have been 
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discussed recently by many authors for example [15]-[19]. The precise measurement of #13 also has impact 
on the global analysis of the mixing parameters. Many groups have already performed global analyses 
recently [8] [9]. We pay particular attention to the result of [8]. There is an interesting indication that the 
mixing angle $23 deviates from the maximal mixing(i.e.,#23 < f (at < 3cr in NH and < 2o~ in IH). And a 
weak hint at the CP-violation phase 8 ~ n. But a reliable result on 8 is not available without the future 
long-baseline neutrino oscillation experiments. 

In the phenomenological view, the new global analysis result may change the form of the PMNS matrix 
significantly. If the value of 623 is much smaller, the absolute value of U T \ (\U T \\ = | ^ ! 2 ^23 _ Ci2 c 23>?i3e i5 l) 
will also decrease. This makes the Hermitian (or symmetric) ansatz of PMNS matrix seem more possible. 
In fact, a Hermitian(and symmetric) ansatz has been discussed by some authors [10][11]. The phenomeno- 
logical view of the Hermitian PMNS matrix also appears in a recent paper [12]. Encouraged by the global 
analysis result [8], we study the possible view that the PMNS matrix is Hermitian in the leading order. 
Consider the strategy that the real PMNS matrix has the following structure [6], 

Upmns =(U,, + AU)K, (6) 
where Uh is the Hermitian ansatz in the leading order, and AU is the perturbation. 

We propose a Hermitian anzatz according to the recent global analysis data in section 2. The relation 
between the possible perturbation and the mass matrix of lepton is discussed in section 3. And we calculate 
next-to-the-leading and next-to-next-to-the-leading corrections. We find that it has a prediction consistent 
with the data very well. Further discussion of the neutrino mass matrix and conclusion appears in section 4. 

O. THE ANSATZ 

To get a proper ansatz, we turn to the experiment data [8]. We simply use the global analysis best-fit 
data in the standard parametrization (3). As observed by Yoni-Ben Tov and Zee [12], the exact Hermitian 
mixing matrix would predict the conservation of CP. So we expect the CP-violation is very small if the 
Hermitian ansatz is proper. Though there is a hint that 8 ~ n, the Hermitian ansatz seems more probable 
when the 8 ~ 0. If the CP- violation originates from a small perturbation of an ansatz, it is more natural to 
assume 8 ~ 0. So we just ignore the possibility that 8 ~ 7r,and use the NH data [8] 

sin 2 6>i2 - 0.307, sin 2 6» 13 = 0.0241, sin 2 6> 23 = 0.386, 5-0 (7) 
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Then the matrix U exp is 



U exp ~ 



0.822 0.547 0.155 
-0.514 0.599 0.614 
0.243 -0.584 0.777 



(8) 



The minus sign of U exp n\ and U ex „-fi can be changed by making a basis transformation 



and 



e L -> diagi l, l, e in W, 
v L -> diag | e i7r , 1 , 1 j v L 



(9) 



(10) 



a phase of K would change, if the neutrinos are Majorana. Then the matrix U exp would look proximately 
Hermitian 



-0.822 0.547 0.155 

0.514 0.599 0.614 (11) 
0.243 0.584 -0.777 

Encouraged by the data, we put forward an simple ansatz of the matrix U. The form of this matrix is the 
following: 



( \ 
-a b c 



U 



(0) 



(12) 



b b b 
c b -a 

This form of the mixing matrix also exhibits symmetry in the lepton sector. The general real and Hermitian 
matrix is 



-a b c 



U 



(0) 



(13) 



b d e 
c e f 

If the Z2 transformation v\ <-> vs and e <-> t keeps the mixing matrix invariant, we have -a = f and b = e. 
We just need d = b to get the form (12), it turns out that this requirement results from a certain form of the 
mass matrix, as will be discussed in sect.4. The unitary property of matrix U helps us to fix the parameters 
a, b and c. The numerical result is 



U 



(0) 



1 



1 



1 



V3-3 V3 V3+3 

1 1 1 

V3 V3 V3 

1 1 1 



V3+3 V3 V3-3 



(14) 



uF has a simple and beautiful structure. Comparing with the parametrization (3), we get #13 * 12.2° 
deviating from experiment result about 3° in this simple mixing pattern, and the degeneration of the 623 and 
On, and of course no CP-violation. According to (6), the deviation of our ansatz from the best-fit data is 



AU 



-0.0333 -0.0304 -0.0563 
-0.0634 0.0216 0.0366 
0.0317 0.00660 0.0117 



(15) 



The magnitude of each element of AU is of 0(0.01 }{U eX p is of O(0. 1)). This implies that all the mixing 
angles should receive the same order correction in a natural way [6]. The relative deviation of 13 is much 
larger than others because of the smallness of #13. We will show in section 3 that it is easy to make a 
perturbation to match the experimental data. 



III. HIGH ORDER CORRECTIONS 

The lepton flavor mixing matrix U pmns is directly related with the diagonalization of the charged lepton 
mass matrix M/ and neutrino mass matrix M v . We take the charged lepton mass matrix to be Hermitian and 
consider the case that the neutrinos are Majorana (There is only the left-hand Majorana mass term). So the 
mass matrix can be diagonalized as 

m e , m„, m T I (16) 

and 

VjM v V v = diagimx, m 2 , m 3 j (17) 

The PMNS matrix is Upmns - Vjv v . In general, it is the product of vj and V v that determines the 
mixing pattern. But the possibility that the leading effect of mixing pattern is caused by V v or V\ with 
the other one as the next-leading order perturbation has been suggested and studied [13]. The struc- 
ture of mass matrix is directly related with the transformation matrix. The mass spectrum of charged 
lepton and neutrino would give some information about the mass matrix structure. The much larger 
mass hierarchy of charged leptons may have some relation with the mass matrix elements hierarchy 
structure. And the mass hierarchy of neutrino is not so large, so it may be natural to think that the 
mass matrix elements are not in the hierarchy structure. One possible and simple charged lepton mass 
matrix structure is that it is nearly diagonalized. Then in the leading order the transformation matrix of 
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charged lepton is close to 1 . So it is the transformation matrix V v that determines the leading mixing pattern. 



We start from the charged lepton mass matrix 





m\ 


e 




Mi - 


e* 




K 




* 

{ T] 


A 


m 3 



(18) 



Where m\,m2 and are all real, e, rj and k are complex. Considering the mass hierarchy of charged leptons, 
we assume the following relation between these parameters. 



|e|, m\ <sc ni2, m^, 



\k\ <§: |e|,mi and ni\ \e\ 



(19) 
(20) 


























1 





111 3 



This charged leptons mass matrix is nearly diagonalized. The leading order structure is assumed to be the 
following. 

l - A 

(21) 

This predicts the mass of electron be zero. This form is reasonable for the big mass hierarchy between 
e and /u « 200). The untiray transformation matrix Vi is trivial, i.e. vj = 1. So the transformation 
matrix V v is just the leading mixing matrix U, . And we assume the transformation matrix V v stays the 
same in the all order calculation. The next-leading order is just related with the structure of the charged 
lepon mass matrix. 



Before we go on to consider the next to the leading order, let us compare the leading order prediction 
with the data. The leading order predicts that sin 2 6^ - sin 2 ^3 ~ 0.349. ^3 i s i n tne 2c" region and of^ is 
just in the 3cr region, ff^ is even worse, since sin 2 8^ ~ 0.0446, which is out of the 3cr allowed region of 
the data, so the next-leading order calculation is necessary. 



Consider (19), the next-leading order mass matrix of charged lepton is [14] 








6 







e* 









,0 





m 3 t 



(22) 



With £ = \e\e' A . This matrix deviates from the diagonalized matrix by off-diagonal element e, but it is also 
easy to be diagonalized by the transformation matrix 



cos 9 sinft?^ 0^ 



V 



(i) _ 



-sinBe 




cos e o 
o 1 



(23) 



Where tan6> - « 0.07, - m T , ni2 = m M - m e , \e\ = ^m e m M and <f> = A. Using the relation 



U h = VY"Vy = V} 1 " UY' , we get the next-leading order mixing matrix: 



/(Dt, 



/(l)tr,(0) 



U 



(1) 



cos 9 
V3-3 



V3 



cos 
V3 



sin Be' 



i.i 



sin 9e 



V3 V3+3 



V3 



cos 6 _j_ singe ljl cosfl _j_ singe ljl cos g _j_ singe ljl 
V3 V3-3 V3 V3 V3 V3+3 



1 



1 

V3+3 V3 V3-3 

There is only one free parameter A here. And the possible CP-violation phase is also included in the 
next-leading order. Now after rephasing (23) and comparing with the standard parametrization (3) it is 
straightforward to write down the mixing parameters after correction, 

1 , 2 



l 



(24) 



sin 2 ^ 



( 



■r 



m e 

— cos A, 



sin 2 d£ ~ 



V3 + 3 3 + 3 V3 V "V 

2 + 2( V3-i)^cosi 
4+ V3 + 2(V3 - 1) .[^cosA 



(25) 



(26) 



sin 2 ^ 1 ' 



2 -% 


^cos^ 


4+ V3 + 2(V3- 1)^ 


f cosi 



(27) 



tan<5 (1) 



~(V3 + l)^f sini 



(28) 



Here we only keep the order 0{ J—)- If A is in the first and fourth quadrant, the degeneration of G23 and 
#12 disappers with 623 becoming larger and 612 becoming smaller, which is in keep with the experiment 
data. The next-leading order correction to sin 2 (r^ is about -0.02, which also can fit the data much better. 
More precisely, 6]~ can lie in the 2cr range if < A < 0.413 or 5.870 < A < 2n, 9\) can lie in the 3<x 



range if < A < 0.687 or 5.596 < A < 2n (see figure 1). ^3 can t> e i n the 2o~ range in a large area of the 
parameter space of A. fits the data very well in the leading order, when A is small or near 27r, #j 2 would 
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lie in the Icr region. The best-fit point can also be included. There are two possible CP-violation phases, 
one is near zero, another is near n. As we have assumed, 6 should be near zero in the natural way, A is 
constrained in the fourth quadrant. The rephasing-invariant parameter J~ « -jj ^J^sinA ~ 0(0.001). 



As we can see, the next-leading order correction can fit the data. But it is still not very satisfied, let us 
consider the next-to-next-to-the-leading order correction now. The mass matrix of charged lepton would be 
of the following structure 



( \ 
m\ e 



M' 



(2) 



6* ni2 
m 3 ) 

It is also very easy to diagonalize by a similar matrix 



(29) 



vf = 



cos£ sin^V p 
-sin£e~' p cos£ 
1 



(30) 



Where tan < 



i £ = l m ° +m ; m3 - , Mr> nt2 + m{ = m - me an d p - A. Besides the parameter A, a new parameter 
appears, i.e. m\. We define m\ = (C 2 - l)m e (with C > 1), where C is real and positive. Then tan£ - C yj~^- 
As we have assume m\ «; |e| and know |e| is of 0( yjm fl m e )~ O(l0)m e , m\ should be order m e , i.e. C is 
of 0(1). It is straightforward to calculate the next-to-next-to-the-leading order correction as the same way 
above 



sin 2 0® =* (— ! — ) 2 1 —- , l^-C cos A, 

13 V3 + 3 3 + 3 V3 V "V 



(3D 



sin 2 6g a 



2 + 2(V3-l)^ 


^CcosA 


4+V3 + 2(V3-l)^ 


^CcosA 



(32) 



sin 2 ^ 



2 ~U 


^Ccosi 


4+V3 + 2(V3-l)^ 


^CcosA 



(33) 



tan5 (2) 



-(V3 + l) A &sin^ 



!-(!+ V3)^fCcos/ 



(34) 
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The new parameter C is a linear decreasing function of 613, so #13 is sensitive to the value of C. It is 
possible to make sin 2 9^ lie in the Icr region, sin 2 9^ can still be in the Icr region in a large interval of 

? (2) 9 (2) 

C and A. sin f? 23 is an increasing function of C, if A is fixed in the first or fourth quadrant. But sin $,3 
increases very slowly, so a large value of C is needed to get the Icr region data. On the other hand, too large 
C would make sin 2 ^ and sin 2 6>j 2) decrease too much. Though A could be adjusted to reduce sin 2 ^, 
the parameter space of A and C is expected to be narrow. We expect the higher order correction will affect 
sin 2 0L . But it is easy to find C and A to keep sin 2 8$ and sin 2 0^ both in the Icr region and sin 2 in 
the 2cr region of data. The CP- violation phase <r ^ is still near zero with a deviation about 10°. Future 
measurement of 8 will test our theoretical consideration. 



We do not consider higher order corrections here. 



IV. THE NEUTRINO MASS MATRIX 



In this last section, we make some comments on the mass matrix of neutrinos. The hierarchy structure 
of charged lepton mass matrix elements give the nearly diagonalized matrix in the leading order. The 
transformation matrix V v is just the product of leading order mixing matrix t/j: ^ and K, i.e. V v = u{ K . 
Using (16), we can get the mass matrix of neutrinos as the following 



My - Vydiag^ mu m 2 , m 3 j V 7 V 



Taking V v =ufK into (34), we get M v 



M v = 



x+y z y+z 
z x 2 y + z 
y + z 2y + z x-y 



With 



(35) 



(36) 



x = ^( w Ji + «?2 + mi) 



(37) 



V3 



(mi - m 3 ) 



(38) 



z = -(m 2 + — = 

3 V3 + 1 



1 



-mi) 



(39) 



V3- 1 

where m\ = m\e lip ,m 2 = m 2 e lia '. From the oscillation experiment and the cosmology observation, we 
know the mass spectrum is m\ « m 2 ,m\{m 2 ) «; or mj, «; m\{m 2 ). So elements of neutrino mass 



matrix do not have the hierarchy structure as the charged lepton mass matrix. But (35) has an interesting 
property. The sums of every column and every row are all equal, 2;=i,2,3 ^v,M - x + 2y + 2z (A - e,/j., r) 
and 2u=e,u,T M v ,m = x + 2y + 2z (i =1,2,3)- This property also appears in some A4 models [20][21], but 
mass matrix (36) show less symmetry than the A4 models mass matrix. 

The logic can also be reversed. If the mass matrix has the structure as (35), we can diagonalize it and 
get the mixing ansatz (13). We expect some new symmetry or underlying mechanism would give the mass 
matrix structure (18) and (36). Then we can know more about the origin of the mass matrix structure, 
which may be related with the high-scale physics. We also notice the different mass form of charged lepton 
and neutrino. But it is not unnatural that the mass matrix origins of charged leptons and neutrinos may be 
due to different mechanisms. 

We show in this paper that a Hermitian ansatz can be made to fit the experimental data very well. 
The mixing ansatz that we propose in section 2 is just a phenomenological consideration, its theoretical 
mechanism is unclear yet. The future precise measurements of the mixing parameters, especially of the 
CP-violation phase 8 in the LBL experiments, will provide more information on the mixing pattern and 
constrain all the parameters. We expect 5 be near zero, and then the Hermitian ansatz can be an important 
hint at the underlying physical theory. 
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FIG. 1: sin #13 as a function of A in the next-to-the-leading order 
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FIG. 2: sin 023 as a function of A in the next-to-the-leading order 
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